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Gapped Z2 spin liquids have been proposed as candidates for the ground-state of the S = 1/2 
quantum antiferromagnet on the Kagome lattice. We extend the use of Projected Entangled Pair 
States to construct (on the cylinder) Resonating Valence Bond (RVB) states including both nearest- 
neighbor and next-nearest neighbor singlet bonds. Our ansatz - dubbed "simplex spin liquid" - 
allows for an asymmetry between the two types of triangles (of order 2 — 3% in the energy density 
after optimization) leading to the breaking of inversion symmetry. We show that the topological 
Z2 structure is still preserved and, by considering the presence or the absence of spinon and vison 
lines along an infinite cylinder, we explicitly construct four orthogonal RVB Minimally Entangled 
States. The spinon and vison coherence lengths are extracted from a finite size scaling w.r.t the 
cylinder perimeter of the energy splittings of the four sectors and are found to be of the order of 
the lattice spacing. The entanglement spectrum of a partitioned (infinite) cylinder is found to be 
gapless suggesting the occurrence, on a cylinder with real open boundaries, of gapless edge modes 
formally similar to Luttinger liquid (non-chiral) spin and charge modes. When inversion symmetry 
is spontaneously broken, the RVB spin liquid exhibits an extra Ising degeneracy, which might have 
been observed in recent exact diagonalisation studies. 

PACS numbers: 



Introduction - Gapped spin liquids are intriguing in- 
teracting quantum spin states which do not bear any lo- 
cal order parameter but possess topological order^ . The 
great excitement for such systems is partly due to their 
potential to realize qubits for future topological comput- 
ers^. Frustrated quantum magnets^ like the spin- 1/2 
Quantum Antiferromagnetic Heisenberg (QAH) model 
on the kagome lattice (see Fig. 1(a)) are believed to be 
good model candidates to host topological spin liquids. 
The Anderson's Resonating Valence Bond (RVB) state^, 
Gutz wilier projected BCS wavefunctions^ or the solvable 
dimer-liquid^ provide simple realizations of such states. 
Nevertheless, Lanczos exact diagonalisations (LED) of 
small clusters^ of the kagome QAH model have revealed 
a large amount of low-energy singlets, a strong signal for 
other competing non-magnetic ground-states (GS), like 
the gapless U(l) spin liquid^ or valence bond crystals 
(VBC) spontaneously breaking lattice translation^^^. 
However, recent Density Matrix Renormalization Group 
(DMRG) studies^ have provided more evidence in fa- 
vor of a topological gapped Z2 spin liquid, in the neigh- 
borhood of a U(l)-Dirac spin liquid state^^ nearby in 
energy^^. 

In bisimplex lattices, like kagome or pyrochlore lat- 
tices, consisting of corner-sharing simplices (e.g. tri- 
angles or tetraedra) located on an underlying bipartite 
lattice, a spontaneous symmetry breaking can be ex- 
pected, resulting in a phase where the left-simplices 
differ from the right-simplices (see Fig. 1(a)). For the 
kagome lattice the inversion symmetry is broken in this 
phase. Interestingly, state-of-the-art LED of kagome tori 
with up to 48 sites^'' have suggested that the GS of the 
kagome QAH model might be a gapped spin liquid with 



8-fold GS degeneracy originating from the product of the 
usual 4-fold topological Z2 degeneracy (on the torus) by a 
2-fold Ising degeneracy associated to spontaneous break- 
ing of inversion symmetry. Following this proposal, we 
construct here a simple ansatz for a topological Z2 spin- 
liquid lacking inversion symmetry that might reflect qual- 
itatively the properties of the GS of the QAH model on 
the kagome lattice. For this purpose we use the Projected 
Entangled Pair States (PEPS) formalism^^ which pro- 
vides a natural construction of RVB wavefunctions^^'^^ 
as well as a simple understanding of their symmetries'^, 
topological and entanglement properties'' , and bound- 
ary theories'^'" ''^. 

PEPS construction - We start with the RVB wave- 
function defined on the kagome lattice of Fig. 1(a) as an 
equal weight (and equal sign) summation of all nearest- 
neighbor (NN) singlet (| tl) - I It)) coverings (NN sin- 
glets are all oriented clockwise on all the triangles). Such 
a state can in fact be represented by a I) = 3 PEPS^^''^ 
(up to local unitaries) in terms of two rank-3 tensors, 
(i) ^ on the sites, and (ii) R\^^^v in the center of 
the left and right triangles, where 5 = 0, 1 are qubits 
representing the two Sz = ±1/2 spin components and 
A, /i, G {0, 1, 2} are virtual indices (see Fig. 1(b)). More 
precisely, A2 ^ = 2 = and zero otherwise, and 
^2,2,2 = 1, and i^A,/i,i. = ^x,fi,u otherwise, with ex,f,,u 
the antisymmetric tensor.'^ One can then group 3 sites 
on each unit cell to obtain a rank-5 tensor (the physi- 
cal dimension is now 2^ = 8) connected on an effective 
square lattice (see Fig. 1(c)). Note that one can arbi- 
trarily choose to group the 3 sites on the left or on the 
right triangles. The amplitudes of the RVB state in the 
local 5'2-basis are then obtained by contracting all vir- 



2 



(a) 












A 







t 



(b) 



c d 









AT =4 p=0.65 Q" 
N=6 p=0.65 / _ 


^ b 




,1,1, 


/7 / 

0.49 /' / " 

I/- 

(a) : 

1 1 1 1 1 1 1 1 1 



O DMRG (S.White) 

■ a=p=0 

a=0.49 p=0.65 



(b) 



-0.02 0.02 0.04 



-0.1 



0.1 



-0.39 
-0.4 
-0.41 ^ 
-0.42 ^ 
-0.43 ! 
-0.44 
-0.45 
-0.46 



FIG. 1: (a) Kagome lattice with a vertical cut (dashed line). 
If the inversion center is lacking, right (light shading) and 
left (darker shading) triangles become non-equivalent but all 
sites (blue dots) keep the same environment. (b,c) PEPS 
representation of the RVB state: three S = 1/2 sites are 
grouped together (b) to construct a rank-5 tensor (c) involving 
2^ = 8 physical states and D = 3 virtual states on the four 
out-going bonds. (d,e) Two dimer configurations around a 
"defect triangle" (in blue). Red ellipses represent singlets 
between two sites. The a parameter (see text) controls the 
relative weights between (d) and (e). 



tual indices (except the ones at the boundary of a finite 
system). This RVB state studied in details in Ref. 22 is 
perfectly isotropic i.e. the energy densities are identical 
on all bonds (if the Hamiltonian itself is isotropic). It is 
also equivalent to a projected BCS wavefunction^^ . 

In the case of the NN Heisenberg S=l/2 antiferromag- 
net on the kagome lattice, "defect triangles" of Fig. 1(d) 
are energetically costly. The above RVB wave function 
has a fixed proportion (1/4) of defect triangles character- 
ized byA = /i = i^ = 2on the 3 bonds of the correspond- 
ing i?-tensor. However, as shown by Zeng and Elser^^, 
the variational energy can be drastically lowered by al- 
lowing local fluctuations around each defect triangle, in- 
volving next-nearest neighbor (NNN) singlets as shown in 
Fig. 1(e). Such an improvement can be performed easily 
within the PEPS formalism by acting with the operator 
I — aFs/2 on every left triangle (see Fig. 1(b)), where I is 
the identity operator, P3/2 is the projector on the fully 
symmetric subspace of three spins 1/2 and a is a varia- 



■ 3/2 



can be re-written 



is the local Hamiltonian 



tional parameter. Note that I 
as (1 - f )I - ^H<, where EI< 
of the left triangle, and does not affect triangles with a 
dimer. Further optimization can be performed by mon- 
itoring the respective proportion of "dressed" left and 
(undressed) right defect triangles, using a different R' 
tensor on the right triangles, i?2 22 = l~/^ (0</3<l) 
and R'^ ^A,/i,zy otherwise. 

Energetics and inversion symmetry breaking- The vari- 
ational energy of the modified PEPS RVB wave function 
can be optimized w.r.t. a and P considering infinitely- 



FIG. 2: (a) Variational energy vs energy difference in left 
and right triangles (varying a, at fixed /3 = /3*). The two 
parabolic "branches" are obtained by interchanging the role 
of left and right triangles in the PEPS construction giving 
two orthogonal wavefunctions and |$r). (b) Energy of 

the (fixed) optimized state vs NNN AF coupling J2 [N^ — 6, 
YC12) compared to DMRG data (S. White). 



long YCn (n = 2Ny) cylinders, i.e. horizontal cylinders 
w.r.t. lattice orientation of Fig. 1(a) as in DMRG stud- 
ies^^ Perimeters with = A (YC8), = 6 (YC12) 
and = 8 (YC16) unit cells are considered. One finds 
a = a* 0.49 and p = 0.65 (finite size ef- 

fects are negligible). Although, by construction, left and 
right triangles are non-equivalent, the anisotropy remains 
quite small at the optimized point. The variational en- 
ergy (per site, in units of the coupling J = 1) for fixed 
/3 = /3* shown in Fig. 2(a) as a function of the difference 
Ae = eR — eL between the energy (per site) of the left 
and right triangles, while varying a around a*, reveals a 
parabolic behavior centered around Ae = Ae* ~ 0.01 
corresponding to only ~ 2.4% difference, independent 
of cylinder perimeter. By interchanging the role of the 
left and right triangles in the PEPS construction, one 
can obtain a degenerate branch with opposite anisotropy. 
The minima of these parabolas at ±Ae* therefore de- 
fine a pair of degenerate (orthogonal) RVB states, 
(dressing the left triangles) and |^r) (dressing the right 
triangles), which transform into each other under inver- 
sion symmetry. Note that the variational energy of this 
"Ising doublet", compared in Fig. 2(b) to the DMRG 
data, is significantly lower than the energy of the NN 
RVB state27'28. 

Topological properties - Let us now consider a finite 
thin (horizontal) torus or cylinder, as depicted in Fig. 3, 
of length N}^ and circumference Ny such that N}^ ^ Ny . 
By construction, each component I^l) or \^r) of the 
RVB doublet inherits the same Z2 topological properties 
as the symmetric NN RVB spin liquid^^ or Kitaev's toric 
code^. First, the parity = ±1 of the number of (NN 
or NNN) singlets cut by the vertical (green) closed loop 
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in Fig. 3(a), is a conserved quantity and defines two even 
and odd topological sectors. In addition, one can con- 
sider threading a Z2 flux inside the torus (or the cylin- 
der) along a horizontal loop (dashed line in Fig. 3(a)). 
Four topological sectors can then be deflned by the ab- 
sence {Wv = +1) or presence (W^ = — 1) of a flux and 
by the parity Gy = ±1. On the cylinder, changing the 
parity (flxed by the boundary conditions) is obtained by 
shifting a line of dimers and leaving two extra spinous s 
and s (of opposite spins) at the two ends. This can be 
viewed as inserting a "spinon" line joining its two ends 
(Fig. 3(b)). Similarly, a Z2 flux is added by inserting 
the matrix diag( 1,1,-1) on all the vertical bonds along a 
horizontal line, leaving two "visons" at the cylinder ends. 
Although these four sectors cannot be distinguished lo- 
cally, a spontaneous ss pair (a spontaneous vison pair) 
whose constituents wind around the torus or the cylin- 
der in opposite directions before annihilating, as shown 
in Fig. 3(b), can detect the presence/absence of a Z2 flux 
(measure the parity). Such virtual processes lead to split- 
tings A^v and AE^ between the (variational) energies in 
these four topological sectors as depicted in Fig. 3(a,b). 
Phenomenologically, one expects AE^ and AE^ to be ex- 
ponentially small in the circumference Ny but to increase 
linearly with the torus/cylinder length Nh, i.e. 



AEs = aNhNyex.p{-Ny/ ^spinon), 
AE^ = 67V/,7V^exp(-7V^/^vison), 



(1) 



where ^spinon and ^vison are the spinon and vison tunnel- 
ing lengths. Note that, in contrast to the torus, the GS 
energy in the odd parity {Gy = — 1) sector on the cylin- 
der involves a flnite energy cost 2Eq associated to the 
presence of two spinous at the cylinder ends. Therefore, 
in the limit N^^ Ny 00, at flxed aspect ratio, the RVB 
liquid is only two-fold degenerate on the cylinder while 
it is four-fold degenerate on the torus. Note that for the 
a = 13 = PEPS, a local parent Hamiltonian has been 

found (for which the degeneracy is also equal to 4 on the 
torus). 20^29 

Spinon and vison tunneling lengths in the kagome 
QHAF can be estimated using the improved RVB state 
as a good ansatz. Finite size scalings of the variational 
energies (per site) in the four topological sectors, on the 
same type of YC cylinders as in DMRG and in the limit 
Nh 00 (we set Eq = 0), are shown in Fig. 4(a). Results 
for the splittings shown in Fig. 4(b) as a function of the 
cylinder circumference can be accurately fltted according 
to Eqs. 1. The spinon and vison tunneling lengths are 
found to be quite short, ^spinon — 0.86 and ^vison — 0.67. 
It is interesting that these values are quite close to those 
for a = /3 = 0, ^spinon — 1-0 and ^vison — 0.65, so that 
we believe that further addition of longer range singlet 
bonds will not signiflcantly modify the coherence lengths. 

Entanglement entropy - We now turn to the entangle- 
ment entropy (EE)^^ which provides direct access to the 
total quantum dimension V. We consider a bi-partition 
of the (very long) cylinders in each of the four topolog- 
ical sectors, as shown in Fig. 5(a). The Von Neumann 




FIG. 3: The four topological sectors {Gv = ±1, = ±1) 
of the Z2 spin liquids on a thin torus (a) or on a thin cylin- 
der (b). Zero or one Z2-flux can be threaded through the 
torus/cylinder (dashed green line). In (b) spinous (visons) 
are shown by black (orange) dots. For = — 1 two spinous 
s and s (with opposite spins) are localized at the cylinder ends. 
For each case, a schematic plot of the low-energy spectrum is 
shown. When inversion symmetry is spontaneously broken, 
each level becomes doubly degenerate. 
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FIG. 4: (a) Finite size scaling of the RVB energy in the 
four topological sectors. The energy averaged between even 
{Gv = +1) and odd {Gv = — 1) boundary conditions (dots) 
can be fitted as Coo i G exp {—Nv/^) (dashed lines), (b) The 
splittings in (a) are plotted vs Nv on a logarithmic scale, re- 
vealing two different coherence lengths. 



(VN) EE entropy —Tt{pa In pa}, where pA is the reduced 
density matrix of one semi-cylinder, can be obtained eas- 
ily in the PEPS formalism^^'^^. Results for a = /3 = 
and a 7^ 0, /3 7^ are shown in Figs. 5(b,c). In both 
cases the data can be fitted according to the area law 
SvN{Ny) = cNy — InP with V = 2 with high preci- 
sion, as expected for a Z2 spin liquid^ We note that 
the inclusion of NNN bonds in the improved RVB state 
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FIG. 5: (a) Bi-partitions of the long RVB cylinders into A 
and B semi-cylinders used to compute the EE. By optionally 
inserting spinon and/or vison lines, four disconnected topo- 
logical sectors can be constructed. VN entanglement entropy 
of the symmetric (b) and simplex (c) Z2 RVB liquids in the 
four topological sectors as the function of the cylinder perime- 
ter A^^ . The dashed lines are linear fits based on the averages 
of the four data points for the two largest (YC12 and YC16) 
cylinders providing \nT> with an accuracy of 10~^. 
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FIG. 6: Entanglement spectra for a bipartition of a Nv = 8 
(infinite) cylinder in the Gv = —l,Wv = +1 (a) and Gv = 
+1, Wv — +1 topological sectors, for a — a" and /3 = The 
levels are labelled according to their total spin S (see legends) 
and momentum K along the cut. (c) String representation of 
the edge states of an infinite semi-cylinder. 



leads to a small increase of the coefficient c, i.e. of the 
wave function entanglement. We also observe that the 
difference of the VN EE in the four topological sectors 
vanishes exponentially fast with Ny. These findings es- 
tablish clearly the Z2 topological nature of the NN and 
"simplex" RVB states. 

Entanglement spectrum and edge states - For fractional 
quantum Hall states with bulk gap, it has been conjec- 
tured that there is a deep one-to-one correspondence be- 
tween the entanglement spectrum^^ (ES) of a partitioned 
system, i.e. the spectrum of = — Inpy^, and the true 
edge spectrum of the corresponding A subsystem with 
a boundary. Recently, it has been argued that a simi- 
lar correspondence applies to non-c/iira/ gapped quantum 
spin systems^^. In fact, the PEPS formalism provides the 
proof that Hi) acts on emerging degrees of freedom on the 
edge^^ so that can be viewed as a "boundary Hamil- 
tonian". It is then natural to argue that such (bond) 
degrees of freedom become "physical" when the system 



is "cut" by a real edge and that their effective interac- 
tion is described by the above boundary Hamiltonian. 
We use this conjecture to argue about the existence of 
gapless edge modes on a RVB cylinder. The ES of an 
infinite bi-partitioned cylinder shown in Fig. 6 shows a 
dense accumulation of levels at low (pseudo-) energy. A fi- 
nite size scaling with the cylinder perimeter confirms the 
gapless nature of the ES. Such a feature is robust and 
was also found for the purely NN RVB state^^'^^ The 
(putative) ES-edge correspondence then implies the exis- 
tence of gapless edge modes for an infinite semi-cylinder 
with a vertical boundary or for a finite (long) cylinder 
with two vertical boundaries. Note that, solvable mod- 
els with Z2 topological order, such as the toric code^, 
generically exhibits a gapped ES^^'^^. The existence of 
gapless edge states^^^^ in the toric code in fact requires 
a particular fine-tuned interaction^^. In the RVB topo- 
logical phase, the boundary Hamiltonian associated to 
a vertical edge is an extended t — J model^^ possessing, 
in addition to topological Z2 invariants, translation and 
SU(2) symmetries (^0 representation). One then may 
expect Luttinger liquid-like (c = 1) zero-energy (non- 
chiral) spin and charge modes^^'^^. However, it is not 
clear whether some degree of fine tuning is realized in 
the class of effective t — J model obtained from i^^. Note 
that, as in solvable models, the edge states of the RVB 
liquid should also have a simple string representation^^ 
as schematically sketched in Fig. 6(c). 

In summary, we have extended the previous PEPS con- 
struction of the NN RVB state for the kagome QAH 
model to include fluctuating NNN singlets bonds in the 
vicinity of defect triangles, resulting in a state very close 
to the optimal wave function in the space of translation- 
ally invariant PEPS with bond dimension 3. This novel 
RVB state dubbed "simplex spin liquid" exhibits a weak 
energy difference (~ 2 — 3%) between the two classes of 
triangles resulting in the loss of the center of inversion, 
as suggested by recent LED. In addition to Z2 topologi- 
cal order, Ising (Z2) order is also present in such a state. 
Studying the topological energy splitting of an infinitely 
long cylinder as a function of its perimeter gives access to 
the tunneling lengths of the spinon and vison fundamen- 
tal fractionalized excitations, which turn out to be of the 
order of one lattice spacing. We also argue that gapless- 
energy modes should be present on edges oriented along 
one of the three crystallographic axes. 
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Identying such a Luttinger liquid from its spectrum is te- 
dious due to the existence of both spin and charge zero- 
energy (left and right propagating) modes with different 
velocities. In addition, a precise identification of the Con- 
formal Field Theory levels in the ES spectrum might be 
further complicated by the presence of strong logarithmic 
corrections. 



